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Abstract 

In this talk some of the theoretical aspects of b — > 57 transitions 
are discussed. The focus is on inclusive decays, since these can be 
computed more reliably. Topics covered are (1) perturbative QCD 
corrections, (2) non-perturbative contributions and (3) effects of "new 
physics" in these decays. 



1 Introduction 



Radiative rare B decays have attracted considerable attention in the last 
few years. After the first observation in 1994, by the CLEO collaboration 
I], data have become quite precise so that even a measurement of the 
CP asymmetry in these decays || |J became possible. As far as data are 
concerned, the situation clearly will improve further, after the excellent start 
of both B factories at KEK and at SLAC. 

B — > X s j tests the Standard Model (SM) in a particular way. Since there 
are no tree-level contributions to these processes in the SM, they can occur 
only at the one-loop level. The GIM cancellation, which is present in all the 
FCNC processes, is lifted in this case by the large top-quark mass; if the 
top quark were as light as the b quark, these decays would be too rare to be 
observable. 

Since the SM contribution is small, these decays have a good sensitivity 
to "new physics", e.g. to new (heavy) particles contributing to the loop. 
In fact, already the first CLEO data could constrain some models for "new 
physics" in a stringent way[|TJ. 

The most general effective Hamiltonian describing these decays is given 

by 

#e//=£ C ^> (1) 

i 

where the are local operators 

Oi...q = four-fermion operators 



7 = m b sa^ u (l + -y 5 )bF 
0' 7 = m s sa^(l- l5 )bF^ 
O s = m b sa^T a {l + lb )bG% 

0' 8 ^m s sa^T a (l- 75 )bG% (2) 

and Q are pertubatively calculable coefficients. 

In any new physics analysis of B decays only the coefficients q are tested 
. The decay B — > X s j (and the corresponding exclusive decays) are prac- 
tically determined by the two operators O7 and 0' 7 , and hence these decays 
are mainly testing c 7 and c' 7 . In the SM these two coefficients are 

^-^^ (3) 
V tb V*C 7 m s (4) 



7 32^2tt 2 

where C 7 is a function of (m t /M w ) 2 , which we shall discuss later. 
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Furthermore, the two operators differ by the handedness of the quarks; 
in order to disentangle these two contributions there has to be a handle on 
the polarization of the quarks or of the photon, which is impossile at a B 
factory. Consequently, from b —> 37 alone only the combination \cj\ 2 + \c' 7 \ 2 
can be determined in the near future. 

Once the effective interaction for the quark transition is fixed, one has 
to calculate from this the actual hadronic process. This step is only for the 
inclusive decays under reasonable theoretical control; for exclusive decays, 
form factors are needed, which either need to be modelled or will finally 
come from the lattice. 

For inclusive decays the machinery used is the heavy mass expansionF]. 
Using this framework for the total rate one can establish that (1) the lead- 
ing term as mb —>■ 00 is the free quark decay, (2) there are no subleading 
corrections of order l/m&, (2) the first non-vanising corrections are of order 
1/ml and are given in terms of two parameters. This will be discussed in 
section ||. Additional non-perturbative uncertainties are induced by a cut on 
the photon energy, which is necessary from the exprimental point of view to 
suppress backgrounds. 

Part I deals with the perturbative corrections, part II with the non- 
perturbative ones. In part III, "new physics" in b — > 57 is considered. 

2 PART I: Perturbative Corrections 

The main perturbative corrections are the QCD corrections, which are sub- 
stantial. These corrections are calculated using an effective field-theory 
framework. To set this up, we have to write down first the relevant effective 
Hamiltonian as in ([I]). The operators appearing in ([I]) mix under renormal- 
ization as we evolve down from the My/ mass scale to the relevant scale, 
which is the mass of the b quark. The coflicient functions are calculated at 
the scale /i = My/ as a power series in the strong coupling 

Cl {M w ) = cf\M w ) + ^±cf\M w ) + -.. (5) 

7T 

Changing the scale [i results in a change of the coefficient functions and in the 
matrix elements, such that the matrix element of the effective Hamiltonian 
remains /i-independent. This change can be computed perturbatively for 
sufficiently large /i, using the standard machinery of renormalization group, 
which involves a calculation of the anomalous- dimension matrix that de- 
scribes the mixing of the operators (0). 

1 A non-exhaustive selection of revies is [p[ [7[ p|, p|. pl| . 
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The solution of the renormalizaton group equation yields the coefficient 
functions at some lower scale //, which take the form (schematically) 





(M 2 W 












V 2 J 



c^)=cf , \M w )^[^ln(^)) (6) 

n=0 




where the b n are obtained from the solution of the renormalization group 
equation. 

The last step is to compute the matrix elements of the operators at a scale 
fi « 771&. This can be done for the inclusive case using the l/m& expansion. 
For the exclusive case, one would need the form factor in the corresponding 
approximation, which cannot be done with present theoretical techniques. 

At present, the leading and the subleading terms of the coefficients have 
been calculated 111], [12], [13], KM, including electroweak contributions p5[, the 



main part of which is due to the correct scale setting in a em . 

A complete and up-to-date compilation can be found in JIE| . Without 



going into any more detail we only quote the result from fig 

Br(B -> A s7 ) = (3.29 ±0.33) x HT 4 . (7) 

where this result includes a cut on the photon energy at E 7tmin = 0.05 m^. 

The QCD corrections are in fact dramatic; they increase the rate for 
b — > 57 by about a factor of two. For example, already at the leading- 
log level we have c^(mb) / '^(Mw) = 1.63. Another indication of this fact 
is a substantial dependence of the leading-order result on the choice of the 
renormalization scale /i. This is usually estimated by varying the scale /i 
between and 2mt,. In this way one obtains a variation of 5 M = ^04% 

for the leading-order result. 

Taking into account the subleading terms reduces the scale dependence 
substantially. In fact, one has at subleading order [l^j 5^ = ^32%, which is 



smaller than naively expected fI7 |. It has been argued that the smallness of 
S/j, is accidental However, arguments have been given recently |18| that 



these cancellations are not accidental. In fact, most of the large radiative 
corrections may be assigned to the running of the b quark mass appearing in 
the operator O7. 

3 PART II: Non-Perturbative Corrections 

Non-perturbative corrections arise from different sources. We shall consider 
here 
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• Long-distance effects from intermediate vector mesons B — ► J/^>X S — > 

A s7 , 

• Subleading terms in the heavy mass expansion: 1 / and 1 / m c correc- 
tions, 

• Non-perturbative contributions to the photon spectrum ("shape func- 
tion"). 

3.1 B J/^X S X s7 

One long-distance contribution comes from the process I? — > X S J and 
the subsequent decay of the (off-shell) J/^ into a photon. The first process 
B — > X s J/ty has a branching ratio of order 1%, at least for an on-shell J/^f. 
Assuming that this is similar for the off-shell case, we have to multiply it with 
another factor 1/m 2 for the propagation of the J/^f and a factor fj/y, since 
the J/^f has to annihilate into a photon. This leads us to the conclusion that 
this contribution is indeed negligibly small. However, one has to keep in mind 
that some extrapolation from q 2 = m 2 ^ to q 2 = is involved, assuming that 
this will not lead to a strong enhancement. 

3.2 1/mb and l/m c corrections 

A set of "standard" non-perturbative corrections arises from the heavy mass 
expansion || [?|, |8], |^, [T(J . As far as the total rate is concerned, we have the 
subleading corrections of order l/m|, which are parametrized in terms of the 
kinetic energy Ai and the chromomagnetic moment A2 defined by the matrix 
elements 

2M H X 1 = (H(v)\h v (iD) 2 h v \H(v)) (8) 
6M H X 2 = {H(v)\h v a^iD^D u h v \H(v)) . (9) 

In terms of these two matrix elements the total rate reads at tree level up to 
order 1/m 2 

r ^M"'W(n^,.). do, 

This result is fully integrated over the photon energy spectrum. One can 
also compute the energy spectrum of the photon within the l/m& expansion, 
which is given, again at tree level, by 

dT G 2 F am\ 
dx ~ 32tt 4 



v tx v th *nc 7 



(ii) 
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Figure 1: Interference between O7 and one of the four-fermion operators 
(0 2 ), leading to a contribution of order 1/m 2 . 



6(1 



Ai + 3A 2 v Ai „ 
=—5 {l-x) + —^6 (1 - X) + ■ ■ 



2m( 



which can only be interpreted in terms of moments of the spectrum. We 
shall return to this point in the next subsection. 

If the charm quark is also assumed to be heavy, one may discuss the 
charm- mass dependence in terms of a l/m c expansion fll9fl . The relevant 
contribution originates from the four-fermion operators (e.g. the operator 
O2) involving the charm quark, see fig. |l|. Expanding the matrix element of 
O2 in powers of l/m c we obtain a local operator of the form 
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O l/rnl = s lfl (l- l5 )T a bG: x e^d x F pa 
which can interfere with the leading term O7 (see fig. p. 



12) 



The detailed calculation [20, ^1|] reveals that this contribution is rather 
small 

b^l/ml C 2 A 2 



9C 7 m 2 c 



0.03 



(13) 



3.3 Non-perturbative corrections in the photon energy 
spectrum 

The non-perturbative corrections for the total rate are thus quite small and 
can safely be neglected against the perturbative ones. However, to extract 



6 



the process B — > X s ~f there has to be a lower cut on the photon energy to get 
rid of the uninteresting processes such as ordinary bremsstrahlung. Clearly 
it is desirable to have this cut as high as possible, but this makes the process 
"less inclusive" and hence more sensitive to non-perturbative contributions 
to the photon-energy spectrum. 

Since we are dealing at tree level with a two-body decay, the naive cal- 
culation of the photon spectrum yields a 5 function at partonic level and 
the l/m b corrections are again distributions located at the partonic energy 
E 1 = m b /2, see fllTl). Clearly (|TT|) cannot be used to implement a cut on the 
photon energy spectrum, since this is not a smooth function. 

The perturbative contributions have been calculated and yield a spectrum 
that is mainly determined by the bremsstrahlung of a radiated gluon. This 
part of the calculation is fully perturbative and enters the next-to-leading 
order analysis described in part I of this talk. In particular, the partonic 5 
function smoothens and turns into "plus distributions" of the form 



dT a s 

dx 71 



'ln(l 




x 



(14) 



where the ellipses denote terms that are regular as x — > 1 and contributions 
proportional to 8(1 — x), which are determined by virtual gluons. 

Here we shall focus on the non-perturbative contributions close to the 
endpoint. The general structure of the terms in the l/m b expansion is 



dT 

dx 



J>(— )>(l-x) (15) 
+0((l/m b y +1 5^(l-x)) , 



where <JW is the ith derivative of the S function. 

It has been shown that the terms with 8^(1 — x)/m\ can be re- 

summed into a non-perturbative function such that the photon energy spec- 
trum becomes 

fx = ^^\ V ^\ 2 \ C ^f^ 1 ~ A) , (16) 

where the non-perturbative fuction / is formally defined by the matrix ele- 
ment 

2M B f(u) = (B\Q v 6(uj + iD + )Q v \B) . (17) 

Here D + is the light-cone component of the covariant derivative, acting on 
Q v , which denotes a heavy-quark field in the static approximation. 
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The shape function is in fact a universal function, which appears for any 
heavy-to-light transition in the corresponding kinematical region. In general 
these transitions should be written as a convolution of a (perturbatively 
calculable) Wilson coefficient and the non-perturbative matrix element 

dT = JdojC (u)(B\0 (u)\B) (18) 

with 

OoM = Q v S(u + iD+)Q v (19) 

At tree level this leads to a simple and intuitive formula in which the mass 
m b is replaced by an "effective mass" m\= m b — u> such that 

dT = Jdu dT tree (m b -> m* b )f(u) (20) 

Since this function is universal, it appears in the semileptonic b — > ulv 
transitions as well as in the b — > S7 decays. At leading twist, this leads to 
a model-independent relation between these inclusive decays, which may be 



used to obtain |14 S /V^&| HHh Hi- 



Moments of the shape function can be related to the parameters describ- 
ing the subleading effects in the l/m b expansion. One has 

dcuf(cu) = l, Jduuf{u) = 0, (21) 
duu 2 f(oj) = — — , I duo uj 2 f / J 



3m 2 '7 3m 2 



Radiative corrections can be included using the machinery of effective 
field theory as described in part I. However, here some ambiguity arises from 
the appearance of a double logarithm (see fll^D), which makes the matching 



ambiguous. Various authors |2q , |26| have used the mass convolution formula 
(p0|), although this has not yet been proven to be correct beyond the tree 
level. 

Finally one may also try to resum the subleading terms in l/m b , i.e. the 
terms of order 5^(1 — x)/m % b +1 in fll5|). This has been discussed in |27|, where 
it has been shown that the relevant operators are 

0»{oj) = Q v {iD^5{iD + + u)}Q v (22) 
0%(u) = iQ v [iD^ 5{iD + + uj)} Q v 

Og , {u> u u 2 )= (23) 
QJ(iD + + uj 2 ) {iDl, iDl} 5{iD + + u x )Q v 

Og'iu^uz) = 

g s Q v 5(iD + + u 2 )G^5(iD + + iOx)Q v , 
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Ey(GeV) 

Figure 2: Partially integrated rate normalized to the leading twist result. The 
three lines with a peak correspond to p2 = (500 MeV) 3 , and A = 570 MeV (solid 
line), A = 470 MeV (short-dashed line) A = 370 MeV (long-dashed line). The two 
lines with a dip have p 2 = -(500 MeV) 3 and A = 470 MeV (dashed line), A = 370 
MeV (dotted line). 

plus the corresponding ones where a Pauli spin matrix appears between the 
quark spinors. 

The effect of the subleading terms can be parametrized by four universal 
functions, which appear again in both b — > ulv and b — > s'j. Using a simple 
but realistic model the effects of the subleading terms may be estimated as a 
function of the lower photon energy cut. In fig. |2] we plot the rate integrated 
from a lower cut as a function of this cut for various values of the parameters. 
As expected, the subleading terms at cut values of 2.3 GeV are of order 10% 
and negligibly small below 2.1 GeV. 

4 PART III: "New Physics" in b — ► 57 

In the Standard Model, b — > 37 is a loop-induced process; it thus has consid- 
erable sensitivity to new physics effects. However, as already pointed out in 
the introduction, any B physics experiment tests the coefficients c, appearing 
in the effective Hamiltonian ([I]) and thus all the information on new effects is 
encoded in combinations of the low-energy parameters q, which have to be 
computed in the Standard Model with the best possible accuracy. Compar- 
ing this to B decay data, it will clearly be impossible to find clean evidence 
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for some specific scenario of new physics. 

At present, no significant deviation from the Standard Model has been 
observed in B —>■ X s ^ nor in any other B decay. Given that there are pro- 
cesses that are sensitive to new effects, B physics (and b —>■ 57 in particular) 
can contribute to constrain new physics scenarios. 

Keeping this in mind one may try various scenarios of new physics and 
calculate the effects on b — > S7, i.e. calculate the coefficients of the low energy 
effective Hamiltonian ([!]) in specific scenarios. There is an enormous variety 
of models for physics beyond the Standard Model on the market, and is is 
impossible to cover all these ideas. 

For that reason I shall only consider two examples, which are instruc- 
tive and demonstrate the kind of sensitivity one may expect. In the next 
subsection I shall consider the Type-II two-Higgs doublet model and in sub- 
section I shall discuss a few recent papers on supersymmetry with large 
values of tan j3. 



4.1 Two-Higgs-Doublet Model (Type II) 

One popular and consistent way to extend the Standard Model is to add 
one or more Higgs doublets. This can be done in various ways, but one well 
motivated way is to have two Higgs doublets where one doublet gives the 
mass to the up quarks, the other doublet to the down quarks. 

Out of the eight degrees of freedom of the Higgs sector, three are needed 
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Figure 4: The tan/?-M#+ plane; contours indicate different experimental 



values for B — > X s ^- } figure taken from [29 . 



to give mass to the heavy weak bosons, while the other five become physical 
states. In particular, the spectrum contains a charged Higgs boson, which 
appears in the loops relevant to b — > 57. The first analysis of this decay in 
this type of two-Higgs doublett model has been performed in |28 . 



The parameters of this model are the ratio of the two vacuum expectation 
values (usually expressed as tan/3 = Vi/v 2 ), the mass M H + of the charged 
Higgs boson, all other parameters are irrelevant for our discussion. 



In fig. ^ (taken from |29|) the branching ratio of b — > 57 is plotted as a 
function of the charged Higgs mass, for three different values of the renor- 
malization scale fi. 

In fig. |] (taken from ||29|1 ) we plot contours in the tan /3-Mh+ plane for 
different values of the B — > Xg'j branching ratio. From this figure it becomes 
clear that there is no large effect induced by enlarging tan (3. One may derive 
bounds on the charged Higgs mass independently of tan/3; the current bound 
is M H + > 314 GeV at 95% CL [0. This does not yet include the new CLEO 
result 0, which will move the bound to even higher values. 



4.2 Supersymmetry with large tan/3 

If supersymmetry were an exact symmetry, b — > 37 would vanish, owing 
to the cancellations between particles and sparticles [[H]]. This means that 
b — > 57 tests the breaking of supersymmetry. Clearly many different scenar- 
ios for this symmetry breaking can be invented, having complicated flavour 
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Figure 5: The rate for B — > X s 7 versus tan /3; for the values of the parameters 



see text. Figure taken from [33 . 



structure. 

Again I shall pick an example from a recent analysis |32|, |33|. In these 
papers it has been pointed out that B —>■ X s 7 can indeed be enhanced in 
scenarios whith large tan f3. Working in the MSSM with a flavour-diagonal 
supersymmetry-breaking sector, the relevant parameters are the charged 
Higgs mass M#+, the light stop mass mj , the supersymmetric // param- 
eter, and the parameter A t from the sector of soft-supersymmetry breaking. 

For large tan/3, renormalization group methods may be used to resum 
these terms |33j and one may confront these results with the recent data. In 
fig. ^ (taken from |33|) we plot the rate for B —>■ Xg'j as a function of tan/? 
for [A = ±500 GeV; the values of the other parameters are M H + = 200 GeV, 
mjj = 250 GeV, all other sparticle masses being at 800 GeV. 

A similar plot can be made for negative A t , but for the parameters chosen 
here this scenario is already practically excluded. 

Given such a scenario, one may also scan over some range for the parame- 
ters and identify regions that are still allowed by the experimental constraints. 
In fig. |6| such a scan was performed with m~ t2 < < 1 TeV, m -+ < m~+ < 1 
TeV, \A t \ < 500 GeV, all other sparticle masses being 1 TeV. 

Clearly B — > X s ^y places significant constraints on the parameter space 
of certain supersymmetric scenarios; however, these studies have been per- 
formed with a flavour diagonal supersymmetry-breaking sector. An analysis 
witout this constraint can be found in |34| 
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Figure 6: Allowed range in the M H +-M^ 2 plane for different values of tan/? 



for two values of rrif . Figure taken from [33]. 

5 Conclusion 

The inclusive radiative rare decay B — > X s ^ is under reasonable theoretical 
control; the latest theoretical prediction |18| is slightly higher than (|7]) 



Br th {B -> X sl ) = (3.71 ± 0.30) x 10" 4 , 

where the difference originates from a different value for the ratio m c /mj; 
while [16[] use the ratio of pole masses, in [18 m^ IS /m^ ole is is used as an 



"educated guess" of NNLO corrections. 

The latest (combined) experimental result is [II 

Br exp = {B -> X sl ) = (2.96 ± 0.35) x 10~ 4 , (24) 

which is in agreement with theory within 1.6a. 

The theoretical uncertainty is mainly determined by our ignorance of 
some of the input parameters (quark masses, mixing angles) and to some ex- 
tent also by the uncertainty of higher-order radiative corrections. Improving 
the current theoretical uncertainty will be very difficult with current theo- 
retical tools. 
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